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NEUTRON-INDUCED  REACTIONS 
IN  TISSUE-RESIIDENT  ELEMENTS 


1.  INTRODUCTION 

In  1973,  encouraging  results  for  the  neutron  radiotherapy  treatment  of  cancer  were  reported  from 
the  Hammersmith  Hospital  at  the  Conference  on  Particle  Accelerators  in  Radiation  Therapy  [1].  At  the 
same  conference ,  preliminary  studies  with  the  objective  of  establishing  controlled  clinical  trials  in  the 
United  States  were  also  reported  [2,3].  Such  trials  have  been  undertaken  at  the  Naval  Research 
Laboratory  in  collaboration  with  George  Washington  University,  at  Texas  A&M  University  in  collabora¬ 
tion  with  the  M.  D.  Anderson  Hospital,  at  NASA  Lewis  Research  Center  in  collaboration  with  Case 
Western  Reserve  University,  at  Fermilab  in  collaboration  with  the  University  of  Chicago,  and  at  a 
number  of  centers  in  other  countries. 

These  studies  have  shown  that  the  dose  with  fast  neutrons  has  been  found  to  be  apparently  more 
critical  than  that  with  conventional  megavoltage  radiation,  with  the  consequence  that  there  is  an 
increased  demand  for  high-precision,  high-accuracy  dosimetry  in  mixed  neutron  and  gamma  radiation 
fields.  In  order  to  satisfy  this  demand,  standardization  of  dosimetry  has  been  studied  by  having 
representatives  from  a  number  of  centers  perform  measurements  in  a  common,  prescribed  radiation 
field  [4,5].  The  users  of  homogeneous  tissue-equivalent  ionization  chambers  generally  obtained  results 
in  the  INDl  project  which  were  within  5%  of  the  mean  value  for  tissue  kerma  in  free  air  and  absorbed 
dose  in  a  matrix  phantom  [4]  but:  "...This  variation  is  probably  not  acceptable  since  it  is  generally 
necessary  for  radiotherapy  and  radiobiology  to  know  absorbed  doses  to  within  an  overall  uncertainty  of 
5%,  which  requires  the  basic  dosimetry  to  be  significantly  better  than  this,"  [5].  It  is  noted  that  the 
standard  fields  that  have  been  employed  are  produced  by  the  d (ifn)a  reaction  at  low  energies  which 
results  in  nearly  monoenergetic  neutrons  with  an  energy  not  too  far  removed  from  14  MeV,  This 
energy  region  has  been  and  is  being  intensively  studied  because  of  nuclear  fission,  fusion,  and  weapons 
effects  [6]. 

The  neutron  beams  at  some  centers,  however,  have  spectra  which  lie  largely  at  energies  greater 
than,  say,  15  MeV.  The  previously  mentioned  dosimetric  problem  becomes  compounded  at  these 
higher  energies  because  of  the  lack  of  detailed  neutron  cross-section  data.  Thus,  energy-dependent 
corrections  which  might  be  applied  to  the  basic  dosimetric  measurements  cannot  be  applied.  Similarly, 
the  correction  of  "tissue-equivalent"  kerma  to  "tissue"  kerma,  which  arises  because  of  the  different  com¬ 
positions  of  these  materials  [7],  cannot  be  made  because  of  the  same  lack  of  data. 

The  present  report  summarizes  progress  in  the  calculation  of  kerma  for  tissue-resident  elements 
and  neutron  energies  greater  than  20  MeV.  Section  2  summarizes  the  methodology  and  kerma  calcula¬ 
tions  for  hydrogen.  Section  3  contains  descriptions  of  the  theoretical  models  which  have  been 
employed  in  the  calculations.  The  parameterizations  of  the  models  and  results  for  the  elements  C\  N, 
and  O  are  contained  in  Section  4,  while  Section  5  compares  the  results  of  the  calculations  with  a 
dosimetric  experiment  performed  with  "NE-2I3"  plastic.  This  latter  section  also  contains  calculated 
kerma  factors  for  tissue  and  tissue-equivalent  plastic.  Conclusions  and  suggestions  for  further  work  are 
summarized  in  Section  6. 
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2.  ENERGY  TRANSFER  TO  HYDROGEN 

The  only  charged-particle-producing  reaction  of  significance  for  neutrons  with  energies  greater 
than  20  MeV  incident  upon  hydrogen  is  elastic  scattering;  it  is  expected  that  approximately  50%  of  the 
total  tissue  kerma  is  due  to  this  process  18]. 

Conventional  Legendre  polynomial  expansions  for  the  neutron-proton  elastic  differential  cross- 
section,  fx(0),  are  unavailable  in  the  neutron  energy  range  En  >  30  MeV.  There  have  been,  however, 
extensive  energy-dependent  determinations  of  the  elastic  scattering  matrix  elements  [9.10]  which  are 
linear  combinations  of  P/Mtf),  and  in  practice,  /  ^  5.  The  absolute  value  of  the  square  of 

the  scattering  amplitude  thus  obtained  is  proportional  to  o-(tf),  and  either  recurrence  relations  or 
numerical  integration  may  be  employed  to  obtain  the  kerma  due  to  elastic  scattering 

K^(E)  =  7r  NaA  (A  4*  l)-2/r  J ^  (1  —  cos  0)  cr  (0)  dicos  9) ,  (1) 

where  A  is  the  mass  of  the  charged  recoil  relative  to  that  of  the  neutron,  and  NA  is  the  number  of  tar¬ 
get  nuclei  per  gram.  Six-point  Gaussian  quadrature  has  been  employed  for  the  evaluation  of  Eq.  (1) 
because  the  technique  is  applicable  to  the  elastic  and  inelastic  scattering  from  the  other  tissue-resident 
(TR)  elements.  Table  1  contains  values  [11]  for  the  energy  transfer  to  hydrogen  in  the  energy  region 
30  ^  En  <  60  MeV. 


Table  l  —  Energy  Transfer  from 
Neutron  to  Hydrogen 


£„(MeV) 

30  40  60 

K[radsf  GO9  cm"2)] 

44.53  41.82  37.87 

3.  THEORETICAL  MODELS  FOR  COMPLEX  NUCLEI 
A.  Optical  Model 

All  of  the  theoretical  models  employed  for  the  complex  nuclei  make  use  of  the  optical  model  in 
some  stage  of  the  calculation.  Charged-particle  production  in  the  compound  nucleus  model  requires 
optical  mode]  transmission  coefficients,  while  applications  of  direct  reaction  models  require  optical- 
model  distorted  wave  functions. 


In  the  optical  model,  the  interaction  of  a  free  particle  and  the  nucleus  is  represented  by  a  potential 
well  of  the  form 


-Uop,  -  Vfix)  +  iWf(xD)  +  4 iWDaD\\^f(xD) 


+  *2i  —  4/<xw)  -»W 

r  dr 


In  Eq.  (2)  the  signs  have  been  chosen  so  that  the  well  depths  T,  W%  Wn,  and  Vw  are,  by  convention, 
positive  in  the  usual  case.  In  general,  the  well  depths  are  expected  to  be  energy  dependent  and  to  con¬ 
tain  terms  proportional  to  (N-Z)/A,  where  N,  Z,  and  A  arc  the  neutron,  proton,  and  mass  numbers  of 
the  target  nucleus.  The  shape  factors  are 

fix)  -  (1  +  c')  (3.1) 
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with 

x  -  (r  -  r0Ain)/a,  (3.2) 

and  similarly  for  xD  and  xso.  The  quantity  rt)A  u~  characterizes  the  radius  of  the  potential,  a  its 
diffusivity,  and  r  the  relative  coordinate  of  the  free  paiticle  from  the  nuclear  center.  The  quantity  Kcoul 
represents  the  Coulomb  interaction. 

The  parameters  of  the  model  are  customarily  determined  by  the  fitting  of  experimental  scattering 
data.  The  procedures  employed  in  this  work  are  discussed  in  the  succeeding  sections  de\  oted  to  the 
individual  TR  elements. 

B.  Inelastic  Scattering 

Some  excited  states  in  the  TR  elements  have  large  cross  sections  for  inelastic  scattering  compared 
to  those  for  elastic  scattering.  In  these  cases,  the  contribution  of  the  excited  heavy  recoil  to  the  total 
kerma  should  be  included. 


In  the  distorted-wave  Born  approximation  (DWBA),  with  the  assumption  that  the  projectile- 
target  interaction  leading  to  transitions  arises  from  a  departure  of  Eq.  (2)  from  spherical  symmetry,  the 
differential  cross  section  for  the  scattered  light  particle  may  be  written  in  terms  of  the  center-of-mass 
scattering  angle  as 


d(T 

dSl 


2/  +  1 


fd*rxf\k/mr)  E^(f)  \?(k,.r) 


dr 


(4) 


=  cr(9)  =  S  c t  DW(0), 


where  /  is  the  multipolarity  of  the  transition,  the  *#+  and  KJ  are  optical-model  wave  functions  for  the 
potential  t/,  the  quantity  p  is  a  phase  space  factor,  and  5  measures  the  strength  of  the  transition.  In 
those  cases  where  experimental  angula  distributions  are  available,  the  quantity  S  is  determined  by  nor¬ 
malizing  theory  to  experiment. 

By  means  of  standard  kinematic  formulas  [12]  ,  the  angular  differential  cross  section  of  light  parti¬ 
cles  may  be  converted  to  an  energy  spectrum  of  heavy  recoils. 

C,  Analogue  Inelastic  Scattering  Model 

The  main  source  of  direct  reaction  protons  is  emission  from  low-lying  final  states  of  the  residual 
nucleus,  which  are  isobaric  analogues  of  highly  excited,  particle-unstable,  collective  states  in  the  target 
nucleus.  When  produced,  these  protons  contribute  most  to  the  high  energy  part  of  the  proton  spec¬ 
trum. 


It  has  been  assumed  that  these  reactions  may  be  treated  with  the  use  of  a  collective  interaction 
between  the  free  particle  and  target  nucleus,  given  by  [13] 


MJ 


4NZ 


1/2 


dU] 

dr 


»7 («.*>. 


(5) 


where  d  corresponds  to  the  distance  between  the  centers  of  mass  of  the  photon  and  neutron  fluids  in 
the  Goldhaber-Teller  model,  and  is  related  to  a  dipole  phonon  creation  operator  A  ,fm  by 

<C-  (h/ltinuo)''2  4- 


(6) 
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In  Eq.  (5),  the  energy  of  the  dipole  oscillation  is  E*  =  fiu>  and  p  is  the  reduced  mass  of  the  proton  and 
neutron  fluids  in  units  of  thv  nucleon  mass,  while  U\  is  the  isovector  part  of  the  optical  potential.  The 
relevant  amplitude  is  then  given  by 


S'12 


< 


*/ 


dU  i 
0r 


(7.1) 


and  the  cross  section  by 


o-  (0)  *  S<r dw(9). 


(7.2) 


The  optical-model  wave  functions  are  and  angular  momentum  coupling  has  been  suppressed, 
and 


S 


2mE 


NZ 


(8) 


Equations  (6)  and  (7)  have  been  applied  to  the  analogues  of  electric  dipole  transitions  in  C,  N,  and  0 
in  order  to  obtain  proton  emission. 

Similar  calculations  corresponding  to  the  excitation  of  magnetic  dipole  transitions  have  not  been 
attempted  because  the  version  of  the  code  DWUCK  [14]  implemented  at  NRL  has  not  been  modified 
to  include  the  effect  of  the  tensor  force  between  nucleons.  This  potential  is  known  to  play  a  central 
role  in  comparing  both  the  angular  distribution  and  magnitude  of  theory  and  experiment  for  such  tran¬ 
sitions.  Other  versions  of  the  code  contain  this  effect,  and  such  calculations  should  be  performed.  For 
both  types  of  transition,  there  are  two  particles  in  the  final  states  and  the  contribution  of  the  heavy 
recoil  to  the  total  kerma  should  be  included. 

D.  The  (n,d)*  (n,t ),  (n,3He)  Reactions 

The  total  cross  sections  for  the  production  of  d,  r,  and  3He  tend  to  be  smaller  than  the  production 
cross  sections  of  protons  and  alpha  particles.  Typically, 

<r(n,d)  =  0.5 a(n,p)  (9) 


and 

c T{n,t)  =  <r{n,3 He)  =  0.03<r  (n>p).  (10) 

The  (nj)  and  (n,  3He)  reactions  will  therefore  be  considered  negligible.  However,  although  the  deu- 
teron  cross  section  is  also  rather  small,  its  contribution  together  with  that  of  the  heavy  recoil  is  not 
negligible. 

The  (n,d)  reaction  receives  most  of  its  contribution  by  the  direct  pickup  of  a  proton  from  the  tar¬ 
get  nucleus  and  may  be  treated  theoretically  with  a  modification  of  the  standard  DWBA  [15]  This 
modification  includes  an  approximate  treatment  of  the  virtual  disintegration  of  the  outgoing  deuteron 
and  can  reproduce  the  experimental  data  satisfactorily.  The  virtual  disintegration  ol  the  deuteron  is 
accounted  for  by  replacing  U{h  the  deuteron  optical  potential  at  energy  Ett  in  the  center  ol  mass  system, 
by  the  potential 

Ud  =  (!/£>„)  J*  [tylr  +  s/2 ))  +  Un(\t  -  s/2|)],w ^(s).  (11a) 
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where 

D0-  JVs*(s)  V„„(s). 


(12) 


In  Eqs.  (11)  and  (12),  the  neutron  and  proton  optical  potentials  Up  and  Un  are  to  be  evaluated  at  one- 
half  of  the  deuteron  energy;  the  deuteron  wave  function  is  <t>(s ),  and  the  neutron-proton  interaction  is 
Vnp(s).  If  Vnp  is  of  zero  range,  Eq.  (11a)  reduces  to 


Vd-  \up(r)  +  Un(r)\E_E4, 2 


(lib) 


where  the  well  depths,  energy  dependencies,  and  geometries  for  the  proton  term  may  be  obtained  from 
the  neutron  potentials  and  reaction  Q-values. 

The  calculations  have,  in  fact,  been  carried  out  in  a  finite-range,  nonlocal,  volume-concerning 
approximation,  and  the  DWBA  cross  section  thus  obtained,  cr^,  has  been  normalized  by  a  strength 
C2S  which  has  been  obtained  from  independent  nuclear  structure  calculations  [16]. 


E.  Compound  Nucleus  Reactions 

In  the  present  context,  any  nuclear  reaction  which  is  neither  shape  elastic  scattering  given  by  the 
optical  model  nor  a  single-step  direct  reaction  given  by  the  DWBA  is  subsumed  into  the  category  of 
compound  nucleus  reaction.  This  accounts  for  most  of  the  protons  and  alpha  particles  that  are  pro¬ 
duced;  in  particular,  the  12C(tf,n'3a)  reaction  belongs  to  this  class  of  reaction. 

The  model  of  compound  nucleus  reactions  which  has  been  adapted  [17]  is  a  multistep  Hauser- 
Feshbach  model  [18]  with  pre-equilibrium  effects  treated  by  the  master  equation  formulation  of  the 
exciton  model  [19]. 


With  the  indication  of  all  relevant  quantum  numbers  except  energy  by  .  and  of  the 

energy  corresponding  to  state  a  by  ea,  the  cross  section  for  a  nuclear  reaction  leading  from  channel  a 
to  channel  /3  is  schematically  given  by 


<r(a  — •  p)  ~  |7rXa2ra  (€a)|  [r^)  (eg) j  £  Tp  Up)  pp  (e^) 


(13) 


In  this  equation,  K(t  is  the  neutron  reduced-Compton  wavelength,  and  TQ(ea)  is  the  model  "sticking 
probability";  therefore,  the  first  factor  is  the  probability  of  forming  a  compound  nucleus.  Similarly,  the 
quantity  within  brackets  on  the  second  line  is  the  probability  of  decay  of  the  compound  system  with 
emission  of  particle  /3  at  energy  and  the  sum  is  extended  over  ail  states  into  which  the  compound 
system  may  decay.  The  level  density  of  the  residual  nuclear  system  is  allowed  to  deviate  from  a  stan¬ 
dard  equilibrium  form  [20]  by  a  fraction  obtained  from  the  master  equations  that  reflect  the  approach  to 
equilibrium.  The  level  density  is  taken  to  be 

PpUp)  *  [l  -  /(F)]p0O(€e)  +  £/(pA£)  H  (pAtp),  (14) 


where 

/(F)-  YJ(pAE)  (15) 

p 


is  obtained  from  the  master  equations,  the  equilibrium  level  density  is  pp„,  and  II  (p,htcfi)  is  the  den¬ 
sity  of  levels  with  p-particles  and  //-holes  at  excitation  €p. 


5 


IIERLING.  BASSHL.  ADAMS.  AND  FRASER 


4.  RESULTS  FOR  INDIVIDUAL  ELEMENTS 
A.  Carbon 

Because  neutron  elastic  scattering  data  have  not  been  obtained  in  the  course  of  this  work,  an  opti¬ 
cal  potential  was  determined  by  simultaneously  fitting  neutron  and  proton  elastic  and  inelastic  scattering 
data  which  were  available  in  the  literature  [21-24].  For  a  fixed  geometry,  and  spin-orbit  well  depth  Vso, 
the  individual  best-fit  well-depth  parameters  were  then  fitted  to  a  linear  function  of  the  bombarding 
energy  yielding  for  l2C  in  the  energy  region  20  ^  <  60  MeV 

V  —  61.0  -  0.375  £  MeV, 

W-  4.0  -  0.033  £  MeV, 

WD  -  1.25  +  0.0167  £  MeV, 

(16) 

R0  -  Rso  -  1.069/4 Ui  fm,  RD  -  1.354/4 1/3  fm, 

a  -  aso  -  0.647  fm,  aD  -  0.731  fm, 

where  all  quantities  are  defined  in  connection  with  Eqs.  (2)  and  (3).  The  angular  distribution  cr(0) 
predicted  by  the  optical  potential  may  be  integrated  according  to  Eq.  (1)  to  obtain  the  elastic  contribu¬ 
tion  to  the  kerma.  Neutron  elastic  scattering  from  12C  at  a  bombarding  energy  of  18.25  MeV  is  com¬ 
pared  with  the  best-fit  optical-model  prediction  and  that  of  Eq.  (16)  in  Fig.  1;  the  kerma  due  to  elastic 
scattering  calculated  from  Eq.  (16)  is  shown  in  Fig.  2,  where  it  is  seen  that  the  dashed  extrapolation 
crosses  the  experimental  point  [25].  The  effect  of  inelastic  scattering  to  the  4.44  MeV  excited  2+  state 
has  been  included  in  the  calculation  of  Table  2. 


9  (°) 


Fig  1  -  Neutron  elastic  scattering  from  ,2C  The  best  fit 
optical  model  prediction  is  indicated  by  the  solid  curve, 
while  the  dashed  curve  is  the  cross  section  predicted  by 
Kgs  The  experimental  data  are  from  Ref  21. 
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Fig.  2  —  Elastic  contribution  to  the  carbon  kerma  calculated  with  the  cross  section 
predicted  by  Eqs.  06).  The  experimental  point  is  from  the  Legendre  polynomial  fit 
of  Ref.  25,  and  the  dashed  part  of  the  curve  is  extrapolated  from  the  solid  part. 


Table  2  —  Energy  Transfer  from 
Neutrons  to  l2C  Products 
in  Units  of  Rads/(109  cm"2) 


27.4 

39.7 

60.7 

Products^v^ 

1119 

1.51 

H 

a 

0.99 

_ 

. 

_ 

2.31 

2.34 

1.71 

12C 

Total 

3.88 

4.64 

4.15 

The  contribution  of  protons  to  the  total  kerma  has  both  a  direct  part  and  a  compound  nucleus 
part.  The  direct  part  includes  only  the  analogue  of  the  inelastic  scattering  to  the  electric  dipole  giant 
resonance  (GDR)  in  12C  using  the  optical  potential  of  Eqs.  (16),  and  a  final-stale  modification  for  the 
isospin  dependence  of  the  optical  potential  [26],  The  calculations  have  been  normalized  according  to 
Eq.  (8)  with  £*=  18.2  MeV,  a  mean  value  for  the  excitation  energy  of  the  resonance.  Although  one 
could,  in  principle,  treat  the  excitation  of  the  ground  state  of  UB  within  a  similar  framework,  the  direct 
reaction  code  employed  does  not  include  the  effect  of  the  tensor  force,  which  implies  an  uncertainty  in 
the  normalization.  Consequently,  it  has  been  omitted. 

The  compound  nucleus  contributions  of  protons  and  alpha  particles  have  been  treated  according  to 
the  discussion  of  Section  3E,  with  the  inclusion  of  the  (w,«),  (n,p).  (n,a)  An,n  ) t  in,  2«), 
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(w,rt'a),  in,n 'ay),  (rt,/i'a«"),  (n,n'ap),  and  ( n,n'2a )  reactions;  of  these,  the  last,  of  course,  implies  an 
(n,n'3a)  contribution.  This  chain  is  expected  to  be  dominant  for  alpha-particle  emission  [27). 

The  optical  potentials  were  either  obtained  in  the  present  work  or  chosen  from  a  recent  compen¬ 
dium  [28].  The  general  level  density  of  Ref.  20  was  parameterized  according  to  the  systematic  formu¬ 
lae  of  that  work  but  with  the  pairing  energy  chosen  to  be  12/ A1  2  in  all  cases,  and  shell  corrections 
determined  from  a  liquid  droplet  mass  formula  [29]. 

Figures  3,  4,  and  5  compare  the  theoretical  calculations  with  experimental  data  obtained  at  the 
University  of  California,  Davis  (UCD)*  at  a  neutron  energy  of  39.7  MeV.  Because  kinematic  shifts 
constrain  the  angles  at  which  experimental  data  can  be  obtained,  data  weie  obtained  at  only  a  selected 
set  of  angles;  therefore  it  was  assumed  that  the  experimental  spectrum  observed  at  the  largest  angle  was 
entirely  due  to  compound  nucleus  formation  and  was  isotropic  in  the  laboratory  coordinate  system. 
This  spectrum,  given  as  id2a/ dEd n)„  _  „max,  was  multiplied  by  the  appropriate  solid  angle,  and  the 
result  added  to  the  spectrum  obtained  by  angle-integrating  the  observed  data.  It  should  also  be  noted 
that  the  theoretical  calculations  have  not  been  convoluted  with  the  1-MeV  (fwhm)  neutron  beam 
energy  spread. 
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Fig.  3  —  Comparison  of  experimental  and  theoretical 
proton  spectra  from  l2C  at  a  bombarding  energy  of 
39.7  MeV.  The  experimental  data  in  the  form 
id^r/dEdfl)  obtained  from  UCD  were  treated  as 
described  in  the  text. 


•Note  added  in  proof  The  experiment  data  contained  in  this  report  are  preliminary.  Final  data  arc  compared  with  the  theoreti¬ 
cal  curves  in  F  P  Brady  and  J.L  Romero,  A vutron  Induced  Reactions  Tissue  Resident  Elements ,  C  rocker  Nuclear  l  aboratory  and 
Department  of  Physics,  University  of  C  alifornia,  Davis.  The  agreement  between  the  theoretical  curves  and  experimental  data 
there  is  generally  at  least  as  good  as  that  shown  here 
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f  ig.  4  —  Comparison  of’  theoretical  ami  experimental  angular 
distributions  for  the  ground  stale  l2C(w.</)  l,B  reaction  at  a 
bombarding  energy  of  39 7  MeV.  The  dashed  curve  was 
obtained  with  the  use  of  Hq.  H  Ibb  the  solid  curve  was  obtained 
with  the  use  of  l\q.  Ulal  supplemented  by  volume  conservation 
and  nonlocality  corrections 


The  excess  of  calculated  low  energy  protons  compared  with  observation  is  believed  to  be  an 
experimental  artifact.  Necessary  gamma-ray  discrimination  and  target  energy  loss  establish  lower  limits 
for  the  energies  at  which  data  may  be  obtained.  The  plotted  points  in  this  energy  region  are  therefore 
lower  limits  to  the  proton  production  spectrum.  The  peak  in  the  proton  spectrum  at  21  MeV 
corresponds  to  the  excitation  of  the  giant  electric  dipole  resonance  in  l2C  and  illustrates  that  a  direct 
reaction  component  must  be  included  if  the  spectrum  is  to  be  fitted.  Similarly,  the  paucity  of  high 
energy  protons  is  believed  due  to  the  failure  to  include  the  analogue  of  the  giant  magnetic  dipole  slate 
in  ,2C.  Although  the  basic  physical  theory  for  its  inclusion  is  known,  it  has  not  been  implemented  for 
the  present  calculations.  It  is  emphasized  that  the  theoretical  calculations  are  absolute;  they  have  not  been 
normalized  to  experiment.  Because  parameters  have  been  obtained  either  from  systematics  or  from 
independent  experiments,  it  would  appear  that  the  chosen  models  are  fundamentally  adequate  for 
describing  the  data. 

Figure  4  shows  a  comparison  of  zero-range/local  (ZR/L)  and  finite-range/non-local  (FR/NL)  cal¬ 
culations  with  experimental  data  obtained  at  UCD.  The  spectroscopic  strength  obtained  from  the  ZR/L 
calculation  is  in  better  agreement  with  that  of  Ref.  16  than  that  obtained  from  the  FR/NL  calculation. 
The  better  agreement  of  the  FR/NL  calculation  with  the  shape  of  the  experimental  data,  however,  sug¬ 
gests  that  it  is  the  more  reliable  calculation,  and  at  higher  bombarding  energies  it  is  known  to  lead  to 
spectroscopic  factors  in  better  agreement  with  theoretical  values. 
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Fig.  5  —  Comparison  of  theoretical  and  experimental  alpha- 
particle  spectra  from  I2C  at  J9.7  MeV  The  experimeniaJ  data 
from  UCD  were  treated  as  described  in  the  text. 


Figure  5  shows  a  corresponding  comparison  of  the  calculated  and  experimental  alpha-particle  spec¬ 
tra.  In  the  present  case  it  is  essential  that  kinematics  be  treated  in  detail.  Correct  kinematics  have  been 
included  only  for  the  (a, a)  reaction.  The  parameterization  of  optical  potentials  for  the  calculation  of 
transmission  coefficients  requires  further  work,  and  is  probably  responsible  for  the  calculated  cross  sec¬ 
tion  being  too  large  at  the  intermediate  energies.  The  calculated  peak  at  25  MeV  is  presumably  due  to 
the  assumption  that  8Be  has  only  a  single  discrete  ground  state.  The  agreement  between  theory  and 
experiment  is  considered  reasonable  for  the  model  chosen. 

Table  2  contains  the  various  contributions  to  the  carbon  kerma  which  have  been  calculated 
theoretically  from  the  present  models.  The  deuteron  contribution  to  the  kerma  includes  all  states  calcu¬ 
lated  in  Ref.  16,  placed  at  the  experimental  binding  energies  and  normalized  by  the  theoretical  spectros¬ 
copic  factors.  The  energy  transfers  to  both  the  deuteron  and  the  recoil  I0B  nucleus  have  been  taken 
into  account,  The  results  may  be  compared  with  other  calculations  [30,31],  At  27.4  MeV,  the  present 
result  is  in  better  agreement  with  that  of  Ref.  31  than  with  that  of  Ref.  30.  At  40  MeV,  the  present 
result  appears  to  be  in  good  agreement  with  that  of  Ref.  30,  but  at  60  MeV  the  present  result  is  sub¬ 
stantially  lower  than  that  of  Ref.  30.  A  complete  comparison  of  the  cascade  model  of  Ref.  30,  the 
present  model,  and  experimental  cross-section  data  would  therefore  be  of  interest. 

B.  Oxygen 

Good  neutron  elastic  scattering  data  were  unavailable  in  the  literature,  and  it  was  necessary  to 
employ  proton  data  in  order  to  obtain  a  neutron  optical  potential.  Although  it  was  possible  to  obtain 
values  of  x2  lower  than  those  which  have  been  accepted  as  yielding  "best  fits,"  they  are  marginally 
smaller,  exhibit  some  deterioration  in  the  visual  fits,  and  possess  somewhat  arbitrary  geometrical 
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parameters.  This  is  illustrated  in  Fig.  6,  which  contains  a  comparison  of  experimental  data  at  49.48 
MeV  for  the  elastic  scattering  of  protons  [32]  from  ,60  with  optical  model  fits.  The  best  fit,  which  is 

shown  by  the  dashed  curve,  was  obtained  by  searching  on  the  three  well  depths  F,  and  Wn,  and  on 

the  four  geometrical  parameters  Roy  a.  Rn,  and  while  the  solid  curve  was  obtained  by  searching  on 
the  well  depths  only  with  the  geometry  constrained  to  be  [33] 

R0  «  Rw  «  1.14 Amfm,  R»  «  (17) 

a  =  aw  ~  0.68  fm ;  a =  0.45  fm. 

Although  the  dashed  curve  is  a  slightly  better  fit  to  the  data  than  the  solid  one,  it  has  been  obtained  at 
the  sacrifice  of  increasing  the  number  of  search  variables  from  four  to  seven,  and  the  parameters  are 
subject  to  the  uncertainties  associated  with  known  ambiguities. 


Fig.  6  —  Comparison  of  proton  elastic  scattering  at  49  48 
MeV  (32)  from  l60  with  optical  model  calculations.  The 
dashed  curve  was  obtained  from  a  seven-parameter  search, 
while  the  solid  curve  was  obtained  from  a  three-parameter 
search. 


With  the  geometry  constrained  to  the  values  given  by  Eqs.  (17),  and  with  Vw  fixed  at  7.0  MeV, 
searches  on  V,  W,  and  WD  were  performed  for  data  at  49.48  [321,  46.1  [34],  30.1  [34],  and  20.7  MeV 
[35].  The  results  of  these  searches  were  then  fitted  by  least  squares  to  second-degree  polynomials  in 
energy  with  the  results  that 

V  -  58.645  -  0.528E  +  0.00237E2, 

W  -  -2.348  +  0.1381  £  4-  (3.52  x  10  -)£2.  (18) 

WD  -  9.107  -  0.1747E  +  (7.98  x  10  4)£2, 

where  all  energies  are  in  MeV.  It  is  emphasized  that  Eqs.  (18)  apply  to  protons  and  only  in  the  energy 
range  20  ^  £  ^  50  MeV.  The  low  energy  behavior  of  both  W  and  Wly  is  unsatisfactory,  and,  in  fact. 
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for  £  <  15  MeV  the  sign  of  W  is  incorrect.  In  order  to  convert  these  results  to  a  neutron  optical 
potential.  Coulomb  contributions  must  be  subtracted  from  the  proton  optical  potential.  The  final  result 
for  neutrons  in  the  energy  region  20  <  E  <  50  MeV  is  the  average  potential  given  by 

V  -  57.38  -  0.528£  +  (2.37  x  10“  *>£2  MeV, 

W  -  -2.348  +  0.138i£  4-  (3.52  x  10“5)£2  MeV, 

WD  «  9.107  -  0.1747£  4-  (7.98  x  10“4)£2  MeV, 

Vso  *  7.0  MeV,  (19) 

Rso-  \A4AUifm .  Rd-  \A0A 1/3  fmt 

a  —  aio  *  0.68  fm,  aD  =*  0.45  fm. 


The  elastic  scattering  of  protons  at  30.1  MeV  from  l60  is  compared  with  experiment  [34)  in  Fig.  7. 
The  solid  curve  was  obtained  by  searching  a  V\  W ,  and  W0  with  the  geometry  fixed  by  Eqs.  (17).  The 
dashed  curve,  which  is  nearly  as  good  a  fit  as  the  solid,  was  obtained  from  Eqs.  (17)  and  (18),  and 
justifies  the  choice  of  Eqs.  (19)  as  the  neutron  optical  potential. 


I  ig  7  —  Comparison  of  proion  clastic 
scattering  at  30. J  MeV  1.141  wnh  optical 
model  calculations  for  The  solid  curve 

was  obtained  from  a  three-parameter  search, 
while  the  dashed  curve  was  obtained  from 
Fq  (19) 


The  direct  contribution  of  protons  was  calculated  assuming  that  this  fraction  of  the  ,hO(n,p) 
reaction  proceeds  by  means  of  excitation  of  the  analogue  of  the  GDR  in  160,  which  occurs  at  a  mean 
excitation  energy  of  17.2  MeV.  The  optical  potential  of  Eqs.  (19)  was  used  for  the  initial  state,  and  the 
isospin  correction  of  Ref.  26  was  applied  to  obtain  the  final  state  optical  potential.  The  direct  popula¬ 
tion  of  the  IftN  ground  state  would  correspond  to  the  analogue  of  a  magnetic  quadrupole  transition  and 
is  presumed  weak. 
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The  compound  nucleus  contributions  to  proton  and  alpha-particle  spectra  were  also  treated  accord¬ 
ing  to  Section  3F  with  the  inclusion  of  the  (n.n),  (n.p).  (n.a),  ( rt.tty ).  ( n.2n ),  in,np)>  in. not). 
{n,nay),  (n,nan),  ( njiap ).  {n.n  2a).  ( n.n  lay) .  in.nlan),  in,n2ap).  in,n  3«  )  reactions,  of  ’hese,  the 
last  implies  the  ( n,n'4a )  reaction. 

Figure  8  compares  the  observed  proton  spectrum  at  a  bombarding  energy  of  60.7  MeV  with  that 
calculated.  The  calculated  low  energy  peak  in  the  cross  sections  is  typical  of  evaporation  processes 
The  calculated  peak  near  45  MeV  corresponds  to  the  centroid  of  the  4.39  and  4.72  MeV,  J"  *  l  states 
of  16N,  but  the  calculated  cross  section  between  25  and  40  MeV  is  too  small.  The  experimental  data 
suggest  the  presence  of  perhaps  three  well-defined  slates,  but  the  density  of  levels  in  lhO  is  so  high  that 
it  is  difficult  to  determine  which  analogue  states  are  receiving  a  direct  reaction  contribution  without 
higher  energy  resolution  data  with  improved  statistics. 


Fig.  8  —  Comparison  of  theoretical  and  experi¬ 
mental  proton  spectra  from  l60  at  a  bombarding 
energy  of  60.7  MeV  The  experimental  data  ob-  _ 

tained  from  UCD  were  treated  as  described  in  • 

the  text.  3 

E 

TO 


f  (MeV) 


Figure  9  contains  a  comparison  of  the  theoretical  and  experimental  alpha-particle  spectra  obtained 
at  60.7  MeV.  Although  there  is  a  gross  failure  to  reproduce  the  observed  low  energy  evaporation  peak, 
there  is  reasonable  agreement  between  theory  and  experiment  for  alpha  particle  energies  greater  than 
15  MeV.  It  must  be  recalled,  however,  that  proper  multibody  kinematics  have  not  been  employed  for 
the  calculations.  Similarly,  not  all  paths,  originating  in  I70  and  containing  open  reaction  channels, 
which  reach  a  given  final  nucleus  have  been  included  in  the  calculations.  Both  effects  may  be  expected 
to  alter  the  theoretical  results. 

Table  3  contains  the  energy  transfer  to  ,60  in  the  energy  region  of  interest.  In  this  case  the 
present  results  at  27.4  MeV  are  in  better  agreement  with  those  of  Ref.  30  than  with  those  of  Ref.  31. 
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At  the  higher  energies,  the  present  results  are  substantially  lower  than  those  of  Ref.  30,  and  just  as  in 
the  case  for  |:C,  a  detailed  comparison  of  the  two  calculational  methods  would  be  of  interest. 


C.  Nitrogen 

There  is  a  lack  of  recent  elastic  scattering  data  for  nucleons  in  the  20  —  60  MeV  energy  range 
from  UN.  However,  an  analysis  of  two  nucleon  transfer  reactions  with  light  nuclei  [36]  has  yielded 
energy-dependent  optical  potentials  for  proton  scattering  from  both  UN  and  l5N.  The  energy  depen¬ 
dencies  of  the  potentials  are  sufficiently  different  so  that  attempting  to  extract  a  symmetry  potential, 
which  would  be  useful,  for  example,  in  calculations  of  the  (n,p)  reaction,  would  be  unreliable.  The 
work  of  Ref.  36  also  considers  14N,  and  applying  Coulomb  corrections  to  the  well  depth  of  the  potential 
UN-Pl  of  that  work  yields  the  n-  UN  optical  potential  for  the  neutron  energy  region  18^ <^50  MeV, 

V  —  59.84  -  0.35£  MeV 

W  -  -5.1  +  0.17£  MeV 

WD  -  6.1  -  0.035£  MeV 

V\0  -  5.65  MeV  (20) 

R0  =  Ml  Axnfm,  1.36 AUifm,  Rso  -  l.0Au*fm, 

a  =  0.644  fm,  aD  *  0.52  fm,  aso  =  0.53  fm, 

with  the  restriction  that  if  the  neutron  energy  predicts  a  negative  value  for  a  well  depth,  it  is  set  equal 
to  zero. 

The  proton  spectrum  is  treated  in  the  manner  previously  described  in  connection  with  the  other 
TR  elements.  The  direct  contribution  has  been  taken  to  be  produced  by  excitation  of  the 
8.06  MeV  J*  -  1"  state  of  14N,  which  is  the  analogue  of  the  6.09  MeV,  Jn  —  1",  T  —  1  state  of  14C. 
The  compound  nucleus  part  of  the  spectrum  includes  protons  produced  in  the  in,p),  (n,py),  (n,pn), 
(n,2p),  ( n,pa ),  (n,pny),  ( n,pnp ),  ( n,pna ),  and  (n,p2nx)  reactions,  where  .v  may  be  y,n ,  p,  a,  ay, an, 
ap,  2«.  This  particular  chain  was  presumed  to  be  the  most  probable  of  several  possibilities  because  of 
the  positive  Q- value  for  the  single  step  (n,p)  reaction,  and  because  the  later  stages  follow  the  emission 
of  a  (virtual)  triton  that  may  be  thought  of  as  being  formed  by  the  incident  neutron  and  a  quasi- 
deuteron  outside  of  a  UC  core. 

Figure  10  compares  the  calculated  and  observed  proton  spectra  at  a  bombarding  energy  of  39.7 
MeV.  The  calculated  cross  section  is  higher  than  that  observed  for  energies  less  than  10  MeV,  but  it 
appears  likely  that  this  is  due  to  the  low  energy  experimental  cutoff.  The  peak  at  31  MeV  is  due  to  the 
6.09-MeV  excited  state  of  l2C,  and  is  somewhat  larger  than  that  observed.  The  calculated  spectrum  is 
lower  than  that  observed  at  the  highest  energies,  but  for  reasons  previously  mentioned,  there  has  not 
been  an  attempt  to  calculate  the  analogues  of  magnetic  dipole  transitions.  In  the  present  case,  such  a 
transition  would  add  to  the  population  of  the  F  0+.  T  **  1  ground  state  of  ,4C. 

Figure  11  compares  the  theoretical  and  experimental  alpha-particle  spectra  at  the  bombarding 
energy  of  39.7  MeV.  The  alpha-particle  optical  potential  for  the  14N  (rt.«)nB  reaction  was  not  taken 
from  neighboring  nuclei  but  is  one  which  has  had  some  success  with  heavier  nuclei  [37].  Reasonable 
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Fig  10  —  Comparison  of  theoretical  and  experimental  proton  spectra  from 
14N  at  a  bombarding  energy  of  39  7  MeV  The  experimental  data  from 
UCD  were  treated  as  described  in  the  text 


14N(n,<*)X 
E  39.7  MeV 

f  EXPT.  (jTe*  1 76  ±  30  mb 
o  CALC.  oyc  139  mb 


t  ig  11  —  Comparison  of  theoretical  and  experi¬ 
mental  alpha-particle  spectra  at  a  bombarding  ener¬ 
gy  of  39  7  MeV.  The  experimental  data  from  UCD 
were  treated  as  described  in  the  text 
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agreement  with  experiment  is  obtained  except  for  the  calculated  peak  at  15  MeV.  Similar  to  the  case  in 
160,  this  is  possibly  due  to  the  artificial  nature  of  employing  continuous  level  density  formulae  for  sys¬ 
tems  as  light  as  those  considered  during  evaporation  from  the  TR  elements.  Similarly,  inspection  of 
cross  section  data  indicates  the  difficulties  of  applying  general  models  to  the  analysis  of  these  reactions. 

The  deuteron  contribution  to  the  kerma  has  been  obtained  by  including  all  of  the  appropriate 
stripping  states  of  Ref.  16,  placed  at  the  experimental  excitation  energies  and  normalized  by  the 
theoretical  spectroscopic  factors.  The  calculation  includes  the  kerma  of  the  recoil  l3C  nucleus. 

Table  4  contains  the  energy  transfer  to  l4N  in  the  energy  region  of  interest. 


Table  4  —  Energy  Transfer  from  Neutrons 
of  l4N  Products  in  Units  of  Rads/(109  cm-2) 

rsZoifciol 


Products^v^ 

27.4 

39.7 

60.7 

P 

0.93 

1.28 

1.73 

d 

0.93 

1.05 

0.45 

t 

- 

- 

- 

3He 

- 

- 

- 

“He 

0.53 

0.76 

1.35 

!4N 

0.38 

0.30 

0.21 

Total 

2.77 

3.39 

3.74 

5.  KERMA  FACTORS 

In  an  interesting  experiment,  time-of-flight  techniques  and  pulse  shape  discrimination  have  been 
employed  to  measure  the  contributions  of  protons  and  alpha  particles  to  the  kerma  produced  by  a  neu¬ 
tron  beam  employed  for  radiotherapy  impinging  upon  NE-213  plastic  [38).  Table  5  compares  E  e 
experimental  results  with  theoretical  results  obtained  in  the  present  report.  The  kerma  factors  for  pro¬ 
tons  and  alpha  particles  are  given  by 

KP-f„KH  +  fcKc<P) 

and 

Ka  =  fcKc{a), 

where  fn  and  fc  are  the  weight  fractions  of  hydrogen  and  carbon  and  the  factors  KH,  Kc(P)t  and 
Kc( a)  are  from  Tables  1  and  2.  The  theoretical  and  experimental  proton  contributions  are  in  fair 
agreement  with  each  other,  but  the  theoretical  alpha-particle  contribution  is  considerably  larger  than 
that  observed.  The  theoretical  value  is,  however,  in  fair  agreement  with  that  of  Ref  31. 

Table  5  —  Theoretical  and  Experimental  Kerma  Factors 
in  NE-213  Plastic* 


Kerma  Factors 


The  experimental  data  from  Ref  .»re  gnen  m  r.ids/f  m 
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The  theoretical  quantities  of  Tables  1-4  may  be  employed  for  the  calculation  of  kerma  factors  for 
tissue  and  tissue-equivalent  plastic  according  to 

KlE)-lLfLKL(E).  (21) 

where  fL  is  the  weight  fraction  of  element  L  in  the  material  of  interest,  and  KL(E)  is  given  by  the  last 
row  of  the  relevant  table  as  a  function  of  neutron  energy  E.  Models  for  muscle  and  tissue-equivalent 
plastic  are  given  in  Table  6,  and  Table  7  contains  the  final  kerma  factors.  The  kerma  factors  for  tissue 
are  consistently  lower  than  those  for  plastic  primarily  because  of  the  enhanced  deuteron  and  alpha- 
particle  production  from  carbon  relative  to  that  from  oxygen.  The  present  calculations  suggest  that  for 
these  models,  the  plastic-to-tissue  dose  conversion  factor  is  of  the  order  of  0.8  in  the  20-60  MeV  neu¬ 
tron  energy  range. 


Table  6  —  Weight  Fractions  for  Muscle 
and  Tissue-Equivalent  Plastic* 


Element 

Muscle 

Plastic 

H 

0.10 

0.10 

C 

0.12 

0.76 

N 

0.04 

0.04 

O 

0.73 

1  0.05 

*Rel  7  Heavier  trace  elements  have 
been  neglected  in  the  present  work, 
and  it  is  assumed  that  only  the  most 
naturally  abundant  isotope  need  be 
considered. 


Table  7  —  Kerma-to-Fluence  Ratios  for  Tissue 
and  Tissue-Equivalent  Plastic* 


£„(MeV) 

Muscle 

Plastic 

27.4 

6.42 

7.67 

39.7 

6.30 

7.95 

60.7 

5.84 

7.17 

*  Rads/ UO9  cm 


The  data  of  Table  7  may  be  conveniently  summarized  by 

K((E)  -  R(E)Kp(E), 

where 

/?<£)- 0,84-  (5.0  x  10 ~*)EMeV 


(22.1) 

(22.2) 


is  understood  to  be  applied  for  20  <  E  <  60  MeV,  and  Kf  and  Kp  are  tissue  and  plastic  kerma,  respec¬ 
tively. 

6,  CONCLUSION 

Obtaining  kerma  factors  for  C,  N,  and  O  for  neutron  energies  greater  than  20  MeV  js  difficult.  In 
order  to  more  firmly  parameterize  possibly  applicable  nuclear  models,  data  should  be  obtained  with 
higher  neutron  energy  resolution  and  improved  statistical  accuracy  in  order  to  resolve  individual  nuclear 
states  which  receive  direct  reaction  contributions.  Whenever  possible,  surrogate  proton  data,  which 
may  be  obtained  more  readily  than  neutron  data,  may  be  employed,  and  isobaric  spin  invariance 
invoked  in  order  to  obtain  corresponding  neutron  information. 
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It  is  likely  that  because  of  intensity  limitations,  future  neutron  experiments  will  necessarily  be 
inclusive,  but  certain  specific  reactions  such  as  uC(n,n' 3a)  or  l60  ( n,n'4a )  should  be  checked  at  one 
energy  as  exclusive  reactions.  Such  kinematically  complete,  multiparticle  coincidence  experiments  are 
difficult  and  costly,  but  would  serve  to  establish  normalizations  and  increase  confidence  in  calculations. 
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